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Abstract 

This paper studies polar sets of anisotropic Gaussian random fields, 
i.e. sets which a Gaussian random field docs not hit almost surely. The 
main assumptions are that the eigenvalues of the covariance matrix are 
bounded from below and that the canonical metric associated with the 
Gaussian random field is dominated by an anisotropic metric. We 
deduce an upper bound for the hitting probabilities and conclude that 
sets with small Hausdorff dimension are polar. Moreover, the results 
allow for a translation of the Gaussian random field by a random field, 
that is independent of the Gaussian random field and whose sample 
functions are of bounded Holder norm. 
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1 Introduction 

Anisotropic Gaussian random fields arise naturally in stochastic partial dif- 
ferential equations, image processing, mathematical finance and other areas. 
Let X = {X(t)\t G I C IR^} be a centered Gaussian random field with val- 
ues in R d , where / is bounded. We will call X an (N,d)- Gaussian random 
field. The distance in the canonical metric associated with the Gaussian ran- 
dom field is yE [||A(s) — X(i)|| 2 ], where || • || denotes the Euclidean metric. 
Polar sets of Gaussian random fields are investigated in 0] under the assump- 
tions that the components are independent copies of the same random field, 
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that the variance is constant and that ^/E [||X(s) — X(t)|| 2 ] < c\\s — t\\@ 
holds with constants c, (3 > 0. The recent works pf and [2:] consider the 
anisotropic metric 



with H 6)0,1]^ and assume y/E [\\X(s) — X(t)\\ 2 ] < cp(s,t). In addition 
they require for the variance only to be bounded from below. In this paper 
the assumptions on the variance and on the independent copies in the com- 
ponents are substituted by the milder assumption that the eigenvalues of the 
covariance matrix are bounded from below. The random fields in the com- 
ponents neither need to be identical distributed nor independent. Hence, we 
require weaker assumptions on the dependency structure of the components 
of the Gaussian random field than 0, B and 0. It follows from an upper 
bound on the hitting probabilities of X that sets with Hausdorff dimension 
smaller than d — YljLi 1 /Hj are P° iar - Our results allow for a translation 
of the Gaussian random field X by a random field, that is independent of 
X and whose sample functions are Lipschitz continuous with respect to the 
metric p. 

As an application we show that an estimator in which calibrates an 
exponential Levy model by option data, is almost surely well-defined. 

2 Main results 

Let X be an (iV, d)-Gaussian random field. Recall that we suppose the index 
set / to be bounded. We will assume the following two conditions. 

Condition 1. There is a constant c > such that for all s,t £ / we have 



Condition 2. There is a constant A > such that for all t E I and for all 
eeE d with ||e|| = 1 we have E[(X$ =1 ejX^t)) 2 ] > A. 

Condition [1] bounds the canonical metric in terms of the anisotropic 
metric p. Condition [2] bounds the eigenvalues of the covariance matrix from 
below. It excludes, for example, cases where X takes values only in some 
vector subspace. 

We will use a uniform modulus of continuity, see (69) in [7|, p. 167]. We 
restate this result in the next inequality. A weaker formulation suffices for 
our purpose and is proved in the Appendix. Let X be an (A r , c?)-Gaussian 
random field, that satisfies Condition [TJ Then there is a version X' of X 
and a constant c > such that almost surely the following inequality holds: 




(1) 



3=1 



^[\\X{s) - X{tW\ < cp(s, t). 



lim sup sup 

40 s,t£i,p(s,t)<e e 




(2) 
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We will always assume that X is a version, which satisfies ([2]). We 
define by Up p (L) :={/:/-)• K d \ \\f(s) - f(t)\\ < Lp(s,t) Vs,t G 1} the 
L-Lipschitz functions with respect to the metric p. In each direction j the 
functions in Lip p (L) are Holder continuous with exponent Hj. We denote 
by Bp(t,r) := {s G IR^I/^s, i) < r} the closed ball of radius r around t. 

Lemma 1. Let X be an (N,d)- Gaussian random field, that satisfies Con- 
ditions U\ and Then for each L > there is a constant C > such that 
for all t £ I, for all r > and for all functions f G Lip p (L) we have 

mf \\X(s)-f(s)\\<r) <Cr d . (3) 

Proof. For all integers n > 1 we define e n := rexp(— 2 n+1 ) and denote by 
N n := Np(B p (t, r)f]I, e n ) the covering number, that is the minimum number 
of p-balls with radii e n and centers in B p (t, r) n I that are needed to cover 
B p (t,r) n I. We have the inclusion B p (t,r) C njlife ~ r 1/H ^tj + r x ' H i}. 
On the other hand each set n^LiI 5 :/' s j + { e n/N) l l H i] can be covered by a 
single ball with radius e n . Hence there is a constant c\ > independent 
of n such that iV n < l\f=i(( 2rN / e n) {1/Hj) + 1) < ciexp(Q2 n+1 ) where 

We denote by {t-™^ G B p (t,r) fl J|l < i < iV n } a set of points such that 
the balls with the centers {tf^} and radii e n cover B p (t,r) n I. We define 

n+l 

r n := /3e„2 2 , 

where j3 > c is some constant to be determined later. For all integers n, k > 1 
and 1 < i < A^fc, we define the following events 



OO 

^{^-/(ftl^r + En}. (4) 



Z=fc 



AW := (J (J 4*> = A^ 1 ) U (J (5) 



fc=li=l i=l 



where the last equality only holds for n > 2. We will show that the prob- 
ability in ([3]) can be dominated by the limit of the probabilities of the sets 

A (n) 

inf \\X(s) - f(s)\\ < r J < lim P(A (n) ). (6) 

seB p (t,r)nI ) n->oo 

For all s G B p (t, r)nl and all n > 1 there exists i n such that p(s, tf 1 ^) < e n . 
By ([2]) we obtain almost surely 

\\X{s) - X(t[ n) )\\ c 
lim sup sup - < — < 1, 

n— >oo sel r n P 
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where the supremum over s is to be understood such that i n varies according 
to s. Let k £ ]c//3, 1[. Especially there is N such that for all n > N we have 

\\X(s) - X(t™)\\ 
sup - — < K. (7) 

sEl r n 

iV+l 

By going over to a possibly greater constant N we ensure that (1— k)c2~ ~ > 
L. On the event 'm£ S £B p (t,r)ni ll^( s ) — — r there exists sq S B p (t,r)f]I 
such that 

oo 

[|*(-o)-/(*o)||<r + £ r <- ( § ) 

Choose zat such that p(so,t[^ ) < e^. Using (J7j), (|HJ) and the Lipschitz 
continuity of / we obtain 



ii*(€ } )-/Oi 

<\\X{t%>) - X(s )\\ + \\X(s ) - f(s )\\ + ||/(. ) - /(O 

oo 

<nr N + r+ ^ r l + L P( S 0,% ) 

/ = AT + 1 

oo oo 

<ktn + r + 77 + (1 — k)c2 2 < r + r/ 

«=Af+l Z=iV 

and (j6]) is established. 

Trivially we have for n > 2 

P(A (n) ) < p^C"- 1 )) +p( J 4W\^("- 1 )) 
and by ([5]) we have 



i=l 



where i' is chosen such that p(t\ n \t[? ^) < e n _i. We note that for n > 2 

P(Af ) \4 n - 1) ) (9) 

(00 00 N 
||X(tf))-/(^)||<r + X:^||X(^- 1 ))-/(^ 1 ))||>r+ £ n 

l=n l=n—l / 

<P - < c 2 r, - Xltjr 1 ))!! > TV-! - Le^i) 
<P (||^(4 n) ) - /(4 n) )|| < c 2 r 5 - ^(i^-^)!! > (/32t - L)e n -ij , 
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where c 2 = 1+/3 2~ exp(-2 m ). We ensure (/32 2 — L) > by choosing 
/3 > L/2. The idea is to rewrite X(t\ n ^) — X(t^ ^) as a sum of two terms, 

one expressed by X(tj n ^) and the other independent of X(t[ n ^). 

Lip p (L) is invariant under orthogonal transformations. By the spectral 
theorem we may choose new coordinates such that the covariance matrix 
at is diagonal. Then the components of X(i- n ^) are independent. By 
assumption (Jj(s) := [Xj(s) 2 ] > 0. We define the standard normal 
random variables 



Note that E[y(4 n) )y(4 n) ) T ] = H hol ds. If 1 [(Xj(s) - Xj(t)) 2 ) > we 
define 



-*/(*) 



Yj(s,t) :-- 



and Y^(s,i) := otherwise. We further define a matrix rj and a random 
vector Z by 



7] :=E 



Y(tf\t^)Y{tf) r 



We observe that \rjjk\ < 1 and hence in the operator norm \\rj\\ < d. The 
random vectors Z(t[ , t^ 1 ^) and Y(Q ) are independent because the co- 
variance matrix is the zero matrix. By the definition of Y(t^) we see 
that Z(t[ n \t^? ^) and X(t^) are independent, too. We want to bound 
P(4 n) \4 n_1) )- If 4 n) = 4 n_1) then P(^ (n) \4 n_1) ) = holds. Thus we 
may assume that p(if \tp 1} ) > 0. © is bounded by 

P (\\X(t[ n) ) - f(t^)\\ < c 2 r, \\Y(t[ n \tt 1] )\\ > ^fw^? ) 

<F f||X(tW) - /(tW)|| < c 2 r, 11^(^,^)11 + 11^(^)11 > 

<F ^(t ) )-/(4 n) )ll < c 2 r, 11^(^,^)11 > ^ n/2 2 ~ 

+F (\\X(t^) - f(t^)\\ < c 2 r, d\\Y{t<T>)\\ > L ) 
■=h + h. 

Each component of Z(t\ , t[? ^) is a weighted sum of at most c£+l standard 
normal random variables with weights in [—1, 1]. Hence the variance of each 
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component is at most (d + l) 2 . In the following cj with I 6 N will denote 
positive constants. By the independence of X(i- n ^) and Z(tf l \tj? ^) we 
have 

-D... (32 n / 2 -L\ 



h =P (HX^) - < c 2 r) P ( HZ^,^- 1 ))!! > 



2c 



2c / 



<C4 r d exp 



/32^ — L ^\ 8c 2 d(d + l) 2 
((32% - L) 2 \ 



8c 2 d(d + l) 2 T 



By the definition of Y(t^) we have with the abbreviation aj = crj{t[ n ^) 

la < / — e V^? ^7 du 

i«-/(tf ) )ii<c 2 r ) ii(^) fe ii>£2^} V2^y ^•••^ 



</ ( — J e *V CT i "/-M/zr 

To sum it up 

' {02? - L) 2 



•(AW) <P(A( n - 1 ))+c 6 r' i iV n 



exp 



16d(d + l) 2 c 2 



k=2 

oo 



16d(d+ l) 2 c 2 



fe=2 



We choose /3 > max(c, L/2) such that fg^m^ws > 2Q. Then the sum 
converges and the lemma follows by ©. □ 

Let H a denote the ce-dimensional Hausdorff measure, see the definition, 
for instance, in 0, p. 129]. Recall that Q = YljLi V^j with Hj as in the 
definition of the metric p. 



6 



Theorem 1. Let X be an (N,d)- Gaussian random field that satisfies Con- 
ditions^ and{^ IfQ<d, then for each L > there is a constant C > such 
that all Borel sets F C M. d and all random fields Y which are independent 
of X and whose sample functions are all in Lip p (L) satisfy 

P (3s G / : X(s) + Y(s) G F) < CH d - Q {F). (10) 

Proof. By Pubini's theorem it suffices to show for all functions / G Lip p (L) 

P (3s G / : X(s) + f(s) eF)< CH d - Q {F). 

We choose some constant 7 > %d-Q{F). By definition of the Hausdorff 
measure there is a set of balls {B(xi, ri) : I = 0, 1, 2, . . . } such that 

00 00 
Fc\jB(x h n) and ^(2 n )^ Q < 7- (11) 
«=o «=o 

For all j we cut the bounded index set / orthogonal to the j-axis with 
distance (ri/N) 1 / 11 ^ between the cuts. Each piece of / can be covered by a 
single ball of radius r\ in the metric p. Hence there is a constant c% > such 
that / can be covered by at most c%r~[® balls. We apply Lemma Q] to these 
balls. By summing up we obtain 

¥(3seI:X(s) + f(s)eB(x l ,r l )) < c 9 rf- Q . (12) 

By (HH) and O we have 

¥(3seI:X(s) + f(s) eF) 
00 

< ^ P (3s G / : X(s) + f(s) e B{x un )) < c lo7 . 
1=0 

We have P (3s G / : X{s) + f(s) G F) < c 10 H d - Q (F), since 7 > H d - Q (F) 
was chosen arbitrarily. □ 



3 Application 

In this section we apply our results to nonpar ametric statistics. Observations 
in the Gaussian white noise model lead to a Gaussian process to which then 
Theorem [1] is applied. We conclude that a certain estimator is almost surely 
well-defined. Existing results on Gaussian random fields have too restrictive 
assumptions for the application considered here. 

The common approach to the estimation of Levy processes is to take 
advantage of the Levy-Khintchine representation. This often involves taking 
a complex logarithm of the empirical characteristic function, see e.g. [1] 
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and 0]. The branch of the logarithm needs to be taken such that the log- 
characteristic function is continuous. This is referred to as the distinguished 
logarithm. To this end the empirical characteristic function has to be distinct 
from zero everywhere which means that zero needs to be a polar set. For 
instance in the set-up of [l[ the following function is estimated: 



V> : R -»■ <D, v ^ - log (1 + iv(l + iv)FO(v)) , 

where T > 0, O is some L X (R) function modeling option prices and T de- 
notes the Fourier transform. The argument of the logarithm is distinct from 
zero for all v G R and the logarithm is taken such that ifi is continuous with 
V'(O) = 0. The function 0{x) is dominated by an exponentially decaying 
function for \x\ —> oo. Especially xO{x) is integrable. To simplify matters 
we suppose that O is continuously observed as in the Gaussian white noise 
model, cf. For e G L 2 (R) we observe 

dO(x) = 0(x)dx + e(x) dW{x) for all x G R. 

Taking the Fourier transform yields 

/oo 
e ivx e(x) dW(x). 
-oo 

Thus the canonical estimator for ip is 

${v) := - log ( 1 + + iu)J r (d(!))(7;)) for all u G R. 



T 

In [l[ a trimmed log- function is used to bound the real part. But this is of 
no importance for the question whether the estimator i/j is well-defined. We 
require the following condition on e. 

Condition 3. There is a p > 1 such that (1 + |x|) p e(x) 2 dx < oo. 

For example, if e G L 2 (R) and e(x) = 0(\x\~ p ) for \x\ — > oo with 
p > 1, then the condition is satisfied. Condition [3] and Lemma H] in the 
Appendix imply the uniform modulus of continuity @ for a version of 
X{v) := e lvx e(x) dW(x). We will assume that X is a version that 
satisfies ([2]). Thus in the definition of ip the argument of the logarithm is 
almost surely continuous. 

Lemma 2. Let e fulfill Condition^ Let O G L 1 (R) suc/i that xO(x) is 
integrable. Then the estimator ip is almost surely well-defined. 

Proof. We have to show that almost surely the argument of the logarithm 
does not hit zero. The process 1 + iv(l + iv)F{dO ){v) equals 1 at v = 0. It 
suffices to consider the process on R\{0}. We rewrite the process as 

iv{l + iv)(— -+TO(v)+ [ e ivx e(x) dW(x) ) . 

\iv{l+iv) ) 
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We define 



/(„):= _i + FO(v) and X(v) := [ e ivx e(x) dW(x). 

iv(l + iv) y_ oc 

On a bounded index set X is an (1,2)-Gaussian random field. We will apply 
Theorem [1] to X, Y = f and F = {0}. It is proved in the Appendix that 
under Condition [3] there is a constant c > such that for all u,t) £ E the 
inequality 

y/E[\\X(u)-X(v)\\ 2 ] <c\u- v \ min W 2 >V . (13) 

holds. This gives reason to the definition p(u,v) := \u — v\ with H = 
min(p/2, 1) G ]l/2, 1]. Thus Condition [1] is satisfied and we have d — Q = 
2 - 1/iT > 0. 

It remains to show that Condition [2] is fulfilled and that / is Lipschitz 
continuous with respect to the metric p. For e = G L 2 (R) we have tp = ip 
and thus ip is well-defined. We will now show that the covariance matrix of 
X(v) is not degenerated if e ^ G L 2 (R) and u^O. Let e£E 2 such that 
e\ + e| = 1. Then there is G [0, 27r] such that e% = simp and e2 = cos</j. 
Consider X as a Revalued stochastic process. The ltd isometry yields 



E[(e 1 X 1 ( y v) + e 2 X 2 (v)) 2 }=E 



4 2' 

(ei cos(ra) + e2 sin(ux))e(x) dVK(x 



oo 



oo 

(ei cos(ux) + e 2 sin(ra)) 2 e(x) 2 dx 



/oo 
(sm(ip + ra)) 2 e(x) 2 dx > 0. 
-oo 



The function 



/oo 
(sin((/9 + ra)) 2 e(x) 2 da; 
-oo 

is continuous by dominated convergence. On ([— V, — 1/V] U [1 / V, V] ) x [0, 27r] 
it takes a minimum Ay > for V > 0. Hence Condition [2] is fulfilled on the 
index set Jy = [-V, -1/V] U [1/V, V]. 

Since xO(x) is integrable we have that TO is Lipschitz continuous on R. 
l/(w(l+w)) is Lipschitz continuous on sets bounded away from zero. Hence 
/ is Lipschitz continuous on iy. Since Jy is bounded it follows that / is 
Lipschitz continuous with respect to the metric p on Jy. 

Thus we may apply Theorem [1] to the index set Iy = [—V,—l/V] U 
[1/V, V}. Since "H d _ Q ({0}) = we obtain P {3v G Iy : X(v) + /(«) = 0) = 
0. Because V > was chosen arbitrarily the Lemma follows. □ 
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Appendix 



Lemma 3. Let X be an (N,d)- Gaussian random field that fulfills Condi- 
tion [Jl Then there is a version X' of X and a constant c > with the 
following property. If e n J, satisfies Y^n=i e n < 00 f or eac/i i? > then 



\\X'(s)-X'L, u 
hmsup sup — < c 



n^oo S! tel,p(s,t)<e n £n ^l og ( e -l; 

/ioWs almost surely. 

Proof. We define the Gaussian random field Y = {Y(t, s)\t, (t + s) G 1} by 
y(t, s) := X(t + s) - X(t). By Condition □ we have 



d((t, s), (t>, s')) := Vn\\Y(t,s)-Y(t',s'W] 

< ci min(p(0, s) + p(Q, s'), p(s, s') + p(t, t')). 

Let e G ]0,l/4[. We denote by D t the diameter of T e := {(t,s) G I x 
B p (0,e)|t + a G J}, where B p (0,e) := {s G WL N \p(s,0) < e}. By the first 
argument of the minimum D e < 2c\e holds. For 5 G ]0, e[ denote by Nd(T e , 6) 
the covering number of T t in the metric d, i.e. the minimum number of d- 
balls with radius 5 that are needed to cover T e . By the second argument of 
the minimum the covering number satisfies Nd(T e ,5) < C2e®/S 2 ®- By the 
previous estimates we obtain with the substitution y = 5 / {c^ 2 ® yfe) 

/ y/logN d (T e ,6)d6 < ^/2Qcl /2Q Ve / " y/logy^dy. 

Jo Jo 

By choosing a possibly larger constant c\ or C2 we may assume 2c\ = c^ 2 ® . 
The integral is solved by 

J V l °gy~ ld y = ~2 — Er^Vloga;- 1 ) + xy/logx' 1 , 

where Erf(y) = -j= e~ t2 dt. Estimating Erf(y) from below yields some 
constant C3 > such that for all x G ]0, l/2[ 



/ Vlogy 1 dy < c?,x^f\og: 
Jo 



Thus for e G ]0,l/4[ 



I C ^lo g N d (T e ,5)dd < c A eJ\oge- l /i 
Jo 
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holds. By Lemma 2.2 in [(j there is a constant K > such that for all 
u > 2Kc4 : e\/\oge~ l l 2 we obtain 

P( sup ||X'(i + S )-X'(i)||>«) <expf--^- 

Using the Borel-Cantelli lemma the statement follows. □ 

Lemma 4. Let e fulfill Condition [3J Then there exists a number c > 
such that for all u,»eR i/ie stochastic process X(v) = e LVX e(x) dW{x) 
satisfies y/E[\\X(u) - X(v)\\ 2 } < c\u - v pn( P /2,i)_ 

Proof. Condition [3] is satisfied for q := min(p, 2) as well. We conclude that: 
E[\\X(u)-X(v)f] 

2" 

(e iux - e ivx )e(x)dW(x) 



=E 



oo 



/oo 
| e i« _ e i"S|2 e ( x )2 dx 
-oo 
/■oo 

< / min(4, (u — v) 2 x 2 )e(x) 2 d3; 



4e(x) 2 dx + / (u — v ) 2 x 2 e(x) 2 dx 

|x|>2|u-u|- 1 J|a;|<2|M-i)|- 1 



< 



J|x|>2|«-i,|-i V 2 F- V I V 



2k_ u |-i\ 2 -'? 



+ / ( p — ) (it - u) 2 x 2 e(x) 2 dx 

J|a:|<2|«-i;|- 1 V Fl / 

/oo 
|x| 9 e(x) 2 dx 
-oo 

This shows the lemma. □ 
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